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We propose generalized portmanteau-type test statistics in the frequency domain 
to test independence between two stationary time series. The test statistics are 
formed analogous to the one in Chen and Deo (2004, Econometric Theory 20, 
382-416), who extended the applicability of portmanteau goodness-of-fit test to 
the long memory case. Under the null hypothesis of independence, the asymptotic 
standard normal distributions of the proposed statistics are derived under fairly 
mild conditions. In particular, each time series is allowed to possess short memory, 
long memory or anti-persistence. A simulation study shows that the tests have 
reasonable size and power properties. 

1 Introduction 

Recently there has been considerable amount of work devoted to testing the inde- 
pendence or non-correlation of two stationary time series, e.g. Bouhaddioui and 
Roy (2006), Duchesne and Roy (2003), Eichler (2006), Hallin and Saidi (2005, 
2007). Consider two univariate time series 

oo oo 

Xit = '^ajUt-j, X2t = '^bjVt-j, t £ Z, (1) 

where {ut} and {vt} are each a sequence of independent and identically distributed 
(iid) random variables. Our goal is to test the null hypothesis that {Xi^} and 
{X2t} are uncorrelated at all lags. If the joint distribution of the two series is 
Gaussian, then non-correlation is equivalent to independence. Given a realiza- 
tion {XktYt=i, k — ^'2' of length n of the processes defined in ([T]), a popular 
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portmanteau-type statistic is based on the sum of the weighted squared cross- 
correlations, i.e. 



n-l 

Gn= K\j/B^)pl,{j), (2) 

where K{-) is a kernel function, Bn is the bandwidth and Puv{j) is the (empirical) 
cross-correlation of two residual series {ut} and {vt} at lag j. Typically the residual 
series are obtained by prewhitening the two time series separately, which can be 
done by fitting a parametric model, such as an ARMA model to each time series 
[Haugh (1976)]. A multivariate extension of Haugh's (1976) idea can be found in 
El Himdi and Roy (1997) and Pham et al. (2003). However, the asymptotic null 
distribution of the test statistic is invalid if the parametric models are misspecified. 
To avoid possible model misspecifications, a long autoregression can be fitted to 
each time series and the consistency of the test is ensured as long as the orders 
of autoregressions grow properly as the sample size increases [Hong (1996)]. See 
Bouhaddioui and Roy (2006) for a multivariate extension of Hong's (1996) test 
statistic. 

All of the work mentioned above seems to exclude long memory time series. In 
the past two decades, a great deal of research has been conducted on long memory 
time series and its importance has been recognized in both applied and theoretical 
time series literature [Doukhan et al. (2003), Robinson (2003), Teyssiere and 
Kirman (2005)]. In this article, we develop tests for independence between two long 
memory time series, which seem lacking in the literature. Our test statistics are 
frequency domain analogues of Hong's (1996) test statistic, which has been recently 
reformulated by Eichler (2006) in the frequency domain. However, the theory in 
the latter paper was developed under stringent conditions. For example, the author 
assumed the existence of all moments and summability of joint cumulants up to 
all orders, which ruled out the interesting long memory case. 

The portmanteau-type statistic for testing the independence between two time 
series bears some resemblance to that for testing the goodness-of-fit for time se- 
ries models. For the latter problem, Chen and Deo (2004a) formulated a gen- 
eralized test statistic in the frequency domain and extended the applicability of 
portmanteau-type test to the long memory case. The development in this paper 
parallels Chen and Deo (2004a) in that for each time series the prewhitening is also 
done in the frequency domain with the unknown spectral density replaced by the 
estimated spectral density. However, our work differs from Chen and Deo (2004a) 
in two important aspects. First, in our feasible test statistics, the spectral densities 
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of Xfct, k = 1,2, can be estimated under parametric assumptions [see Chen and 
Deo (2004a)], or via a nonparametric approach, where no parametric models need 
to be specified. See Sections 12.11 and 1 2.21 for more details. Second, our results are 
applicable to not only short/long memory time series, but also the anti-persistent 
case. Similar to correlation-based independence tests, the tests proposed can't dis- 
criminate between non-correlation and independence for non-Gaussian time series, 
although our asymptotic theory allows for non-Gaussian linear processes. Never- 
theless, they still help the practitioner extract useful information from the data, 
especially when the null hypothesis is rejected. Finally, we note that the tech- 
niques and results developed here are not directly applicable to relax the Gaussian 
assumption in Chen and Deo (2004a), where the proofs require some sharp bounds 
for the products of discrete Fourier transforms. 

Now we introduce some notation. For a vector x = (xi,-- - iXq)' G M^, let 
\x\ = (X]^=i xfY^'^ . For two sequences (a„), (6„), denote by ~ 6„ if anl^n 1 
as n — > oo. Denote by — >£) and — >p convergence in distribution and in probability, 
respectively. The symbols Op(l) and Op(l) signify being bounded in probability 
and convergence to zero in probability respectively. Let A^(0, 1) be the standard 
normal distribution. 

The paper is organized as follows. In Section [21 we introduce our test statistics 
and derive their asymptotic null distributions. The size and power properties of 
our tests are examined in Section [3] through simulations. The technical details are 
relegated to the Appendix. 

2 The Test Statistics and Their Asymptotic Null Dis- 
tributions 

Throughout the paper, we assume X^j, k = 1,2, admit the linear processes of the 
form ([T|), where the innovations ut and vt satisfy the following assumptions. 

Assumption 2.1. The series {ut} and {vt} are each a sequence of iid random 
variables with mean zero and finite fourth moment. Without loss of generality, 
we assume var(Mt) = vav{vt) = 1. Denote by C4(u) = cum(iio, mq, no, no) and 
C4{v) = cum(i;o, vq, vq, vq). 

Let i = i/— 1 be the imaginary unit. For a complex number c let c be its con- 
jugate. For any two processes {Zit,Z2t,t £ Z}, denote by fziZiW and fz2Z2W 
their spectral densities respectively; define their Fourier transforms, periodograms 
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and cross periodogram by 

I " 

wz,{X) = ^= VZfcte'*^, /Zfc(A) = \wz^{X)\'^, and /ziZ2(A) = wz,Wwz^{\) 



for A: = 1,2. Let Aj = 2tt j/n,j = 1,2, ••• ,n, be the Fourier frequencies. To 
motivate our test statistic, we express ([2]) in the frequency domain. Define the 
residual cross-covariance function by 



^ Er=i+i ^t-jVt (i > 0) 

Er=i-i '^tvj+t (j < 0) 



RiivU) = 

and the auto-covariance function for the residual series ut and vt by 



Ruu{j) = n ^ ^ utut-.\j\, Rvv{i)=n ^ ^ ^'t^i-lih 

t=l+|j| i=l+|i| 

Then pui,{i) = Ruv{j)/{Ruu{0)Rm{0)V/^. Let 

where K(-) is a symmetric kernel function with K(0) = 1. Then 

Gn = 27r r \fui,{X)fdX ( r fuu{X)dX ^ fu,{X)dx) • (3) 



Let W{X) = {2tt) ^ X^|/i|<„ -^(^/^n)e be the spectral window function corre- 
sponding to the kernel function K{-). We can write fuv{X), fuu{X) and fv{,{X) into 
the following equivalent forms in the frequency domain (see Equation (3) in Chen 
and Deo (2004a)), 

fz,Z2{X)= W{X-w)Iz,z,{w)dw, Zi,Z2=u,v. (4) 
Jo 

The expressions ^ and (jH) motivate us to propose the following test statistic: 

Tn = y.\fxxX2{^i)?' { — 2^/XiXi(A/) l^fx2X2{^l)\ ) 

1=0 (. 1=0 1=0 J 



where 



j=i V fxhXh {Xj)fxkXk (Xj) 
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Similar to Chen and Deo (2004a), our test statistic T„ is obtained by discretizing 
the integrals in ([3j) and ^ with replaced by IxiX2 1 \J fxiXi /X2X2 ^iid luu 
replaced by IxixJfxiXi (-^XaXa/ZxaXa)- Since fxiXi, fxiXi and /xaXa 
are evaluated only at Fourier frequencies, our test statistic T„ is mean- invariant. 
Our test statistic T„ is infeasible since we do not know the true spectral densities 
/xiXi(A) and fx2X2W practice. In the next two subsections, we introduce two 
ways to estimate the spectral densities, which lead to two feasible test statistics. 

To establish the asymptotic null distribution of T„, we make the following 
assumptions on the kernel function K{-) and the bandwidth 

Assumption 2.2. The kernel function K : M ^ [—1,1] has compact support on 
[— 1, 1]. It is differentiable except at a finite number of points and symmetric with 
A'(0) = 1. 

The assumption that the kernel function has compact support can be relaxed; 
see Chen and Deo (2004a). Here we decide to retain this assumption to avoid more 
technical complications in view of our long and technical proof. It is worth noting 
that several commonly-used kernels in spectral analysis, such as Bartlett, Parzen 
and Tukey kernels, satisfy Assumption 12.21 (see Priestley (1981), p 446-447). 

Assumption 2.3. The bandwidth Bn satisfies (log^ n)/Bn and i?„(log^ n)/n 
0. 

Let A{\) = E^o%e'^^ and B{\) = E^o^i^*^^- The following assumption 
is made regarding the long memory behavior of {X^t}, k = 1,2 and is satisfied 
by two commonly- used long memory time series models: FARIMA (fractional 
autoregressive integrated moving average) models and fractional Gaussian noise 
[Beran (1994)]. 

Assumption 2.4. For k = 1,2, assume fxkXkW ~ |A|~^'^'="Gfc as A ^ 0, where 
dko £ 1/2) and G (0,oo). Further we assume that 

\dA{X)/dX\=0{\A{X)\\X\^^) and \dB{X)/dX\ = 0{\B{X)\\X\-^) 

hold uniformly in A G (0, vr]. 

For A; = 1, 2, the process {Xf^-t} is said to possess long memory if d^o G (Oi 1/2), 
short memory if d^o = and anti-persistence if d^o G (~l/2, 0). Our results cover 
the short memory and anti-persistent cases as well. 
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Theorem 2.1. Suppose that the two processes {ut} and {vt} are independent. Un- 

>D iV(0,l), 



der Assumptions \2. 1\2.4\ we have 

nTn - BnS{K) 



V^Br^diK) 

where s{K) = J^^K'^{x)dx and d{K) = J^^K^{x)dx. 
2.1 Feasible Test Statistic I 

In this subsection, we estimate the spectral densities in a parametric way, so the 
fohowing parametric assumptions are imposed on fx^x^iX), k = 1,2. 

Assumption 2.5. For k = 1,2, let the spectral density of the process {X^t} be 
fk{X;dko), where is the true parameter vector that lies in the interior of the 
compact set 0^ C M'^*, qk G N. Suppose that the estimator dk satisfies \6k —Oko\ = 
Op(n-V2), k = l,2. 

Here we can take the Whittle pseudo-maximum likelihood estimator as dk- The 
root-re asymptotic normality of Whittle estimator has been established by Hannan 
(1973) for d = 0, and by Fox and Taqqu (1986), Dahlhaus (1989) and Giraitis 
and Surgalis (1990) for d G (0, 1/2). See Velasco and Robinson (2000) for the case 
d £ (-1/2,0). 

Let 6 = (^'1,^2)' ^'^d = {01,02)' ■ Then we can replace /xiXi(A) and fx2X2W 
in Tn with /i(A;^i) and /2(A;^2) respectively and get the following feasible esti- 
mator 

A^2 f n„ 27r ""-^ 1 

n ^-^ n ^-^ re ^-^ 

1=0 I 1=0 1=0 ) 

i=i \J fh{Xj;0h)fk{Xj;0k) 

To obtain the asymptotic null distribution of Tn{0), we further make the fol- 
lowing two assumptions; compare Assumptions 6&:7 in Chen and Deo (2004a). 

Assumption 2.6. For k = 1,2, let = (y^fcc dko)' ■ The spectral density /fc(A; ^^o) = 
/fc(A;4o)5fc(A;^fco)' where and gl are even functions on [-7r,7r], /fc(A;4o) ~ 
c(dfco)A~^'^'''° as A — > for some constant c((ifco) > 0, gl{X;Pko) is bounded away 
from zero and differentiable on [— 7r,7r]. Further assume that the q^th component 
of Gfc is contained in the closed interval [— k, k] for some k G (0,0.5). 
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Assumption 2.7. For k = 1,2, 



1. logfk{X;0k), dlog fk{X] Ok)/ dOk^ and log fk{X; Ok) / dOk^dOk, are continu- 
ous at all {\,0k) except A = 0. Further, 



sup sup \Xr''''fk{\]Ok) 
Ae[o,27r]flfceefc 



Ak, for some Ak G (0, oo) 



2. For any (5 > 0, 



sup sup |A| 

Ae[o,27r]0fceefc 



91og/fc(A;^fc) 



and 



sup sup lAI*^ 

AG[0,27r] Bk^Qk 



log fk{X;Ok) 



dOk^dOk^ 



Ak, for some Af^ G (0,oo) 



Ak, for some A^ £ (0, oo). 



(5) 



(6) 



3. There exists a constant C such that 



(7) 



uniformly in A and all 6^/^^ = (^9^^^', and O'^^'^ 



In Chen and Deo (2004a), the same (and equivalent) form of conditions as 
dSD and © were imposed on df^\X;Ok)/dek^ and f^\X;Ok) /ddk^ddk^. The 
condition ([7]) is a variant of the Assumption 7(iii) in their paper. It is not hard 
to verify Assumptions 12.61 and 12.71 for fractional Gaussian noise and FARIMA 
processes. 

Theorem 2.2. Suppose that the two processes {ut} and {u^} are independent. Un- 
der Assumptions \2. 1^2. T[ we have 



nTnjO) - Bns{K) 
yj2Bnd{K) 



>D iV(0,l). 



In forming our feasible test statistic Tn{0), we assume correct parametric spec- 
ifications for the spectral densities of {X^t}, k = 1,2. In practice, model selection 
technique [Beran et al. (1998)] can be used to identify the best parametric model 
for each time series. 
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2.2 Feasible Test Statistic II 



In this subsection, we shall restrict our attention to the fractionally differenced 
autoregressive models of order oo, denoted as FAR(oo,(i) (i.e. the d-th. fractional 
difference follows an infinite order autoregressive process). Our consideration of 
this class of models is motivated by the results in Bhansali et al. (2006), who ob- 
tained the asymptotic normality of estimated coefficients when fitting an FAR(p, d) 
model (i.e. the d-th fractional difference follows an autoregressive model of order 
p) to the observations from an FAR(oo, d) process with p ^ oo as — > oo. 
For k = l,2, let 

oo 

(1 _ 5)d.oXfc, = y^,, Y.''kj^kt = £kt, (8) 

j=0 

where B is the backward shift operator, Skt are iid mean zero random variables 
with al = Ee|j and Ee^^ < oo. Assume a^Q = 1 and Yl'jLo \ (^kj\ < Then X^t 
has the spectral density 



l_2|i JA|-2dtn 



/x,x,(A) = (27r)-V^|l-e 



2 

oo 
j=0 



k = l,2. 



It is easy to see that the processes {X^t}, k = 1,2 defined in ([8]) are a subclass of 
([1]). The following assumption was also made in Bhansali et al. (2006). 

Assumption 2.8. For some e > 0, 

oo 

afc(x) := ^ aljX^ / 0, \x\ < 1 + e, k = 1, 2. 

j=0 

Remark 2.1. As stated in Remark 2.1 of Bhansali et al. (2006), Assumption 12.81 
implies that Cfc(x) = l/ak{x) has the expansion 

oo 

Ck{x) = ^cljX^, \x\ < 1, cIq = 1, k = 1,2 

j=0 

and that max(a°j , c^^ ) < Cr;^, j € N for some constant C and G (0, 1). It 
is satisfied by FARIMA processes with finite autoregressive and moving average 
orders. 

Following Bhansali et al. (2006), for A; = 1,2, we fit an FAR{pf^,d) model to 
the observations {A'/ct,t = I,-- - ,n} generated from the process ([8]). Let 7^ = 
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(dfc, Ofci, Ofc2, • ■ ■ i^^kp^y be the resulting estimates of the parameter vector 7^ = 
(dfc, afci, Ofc2) • • • > CLkpk)' ™ the FAR(pfc, d) model. Note that the spectral density of 
an FAR{pk, d) model is 

-2 



fk{^;ik) 



27r' 



1 



JA|--2d.. 



j=0 



where a^o = 1. 



Let 7 = (71,72)' and 7 = (7'i,72)'- We form our test statistic Tn{^) by replacing 
fXkXkW in T„ with /fc(A;7fc)> i-e- 

\ -1 



1=0 



— yz fxixii^i) — yzf ^2X2(^1 

1=0 1=0 



where 



n-l 



W{\-\,)Ix,x,{X,) 



h,k = 1,2. 



Since T„(7) is free of u^, we set af, = 1, k = 1, 2. 

Let 7^0 = (c^fcO! flfci) ■ • • ; '^kp^y^ k = 1,2. We require certain convergence rate 
for 7;,, which has been obtained in Bhansali et al's Theorem 3.1. 



Assumption 2.9. Assume that |7^ — 7fcol = Op{y^pk/n), k = 1,2. 

Theorem 2.3. Let Assumptions \2.1^2.4\ and \2.8^2.9[ hold. Suppose that the two 
processes {eu} and {£21} are independent. Assume that for k = 1,2, 



Then we have 



X] \(^kl\ = o(.V^/{n log n)), pk^oo, 

and pI log^ n = o(S„), = o(n). 

nT„(7) - BnsiK) 



(9) 
(10) 



N{0,1). 



yj2Bnd{K) 

Let = Zfc logn, where 1^ > — (logr^)"^. Then ([9]) holds under Assumption 
(see Remark l2.ip and (fTOl) reduces to log^n = o{Bn) and (log^ n)i?„ = o(n). 

Remark 2.2. The frequency-domain prewhitening idea involved in r„(7) is similar 
in spirit to Hong's (1996), where the prewhitening is done by fitting a long autore- 
gression in the time domain for each time series. However, Hong's test statistic is 
valid only in the short memory case, while Tn{'y) is applicable to short memory, 
anti-persistent and long memory time series. 
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The theoretical results presented in this section state the asymptotic null dis- 
tributions of the infeasible test statistic T„ and two feasible ones Tn{6) and T„(7). 
It would be desirable to study the asymptotic distributions of these test statistics 
under either fixed or local alternatives. At this point we are unable to obtain any 
results in this direction due to some technical difficulties. We conjecture that for 
fixed alternatives, the rate at which T„ (T„(^), T„(7)) diverges is n/y/B^; compare 
Theorem 4 in Hong (1996) and Theorem 2 in Bouhaddioui and Roy (2006). In 
addition, our proofs heavily rely on the iid assumptions on the innovations {ut} 
and {vt}- It would be interesting to see to what extent these assumptions can be 
relaxed and if the asymptotic distributions of these test statistics are still valid for 
nonlinear models. 

In the next section we examine the size and power properties of our tests 
through Monte Carlo simulations. 

3 Simulation Studies 

We consider the following two models: 

{I - Bf\l - 0.5B)Xu = Ut, {1 - Bf-\l - 0.5B)X2t = Vt, (11) 

and 

{1- B)~^-^Xu = {l + 0.5B)ut, {I- B)-'^-^X2t = {l + 0.5B)vt, (12) 

where {ut} {{vt}) are iid A^(0, 1) random variables. So in pT]) ( (fT2]) ). Xu has 
moderate long memory (antipersistence) and X2t possesses strong long memory 
(antipersistence) . Two series lengths (n = 64, 128) are investigated. 

The simulation of two FARIMA series from model (jlip involves the following 
steps: 

1. Generate two iid innovation sequences {ut}t=-4ooo,-- ,n and {vt}t=-iooo, - ,n 
with pre-specified correlations. 

2. Simulate two AR(1) time series Xu = 0.5Xit_i +ut and X2t = 0.5X2(_i +vt 
recursively with the first 1000 starting values subsequently discarded. 

3. Generate two FARIMA series by applying a truncated autoregression (corre- 
sponding to the fractional differencing) of 3000 lags [see Martin and Wilkins 
(1999)] to {Xit}t=-3000,- ,n and {X2t}t=~3000,- ,n respectively. Finally, we 
only retain the last n observations to ensure good approximation. 
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The simulation based on model (jl2p is similar except that we generate two 
MA(1) time series in the second step. In the calculation of Tn{d), the Whittle 
estimator 9 was obtained assuming an FAR(1, d) model for each time series. Thus 
the assumption is incorrect if the data are generated from model (jl2p . This will 
give us some idea of the consequences of model misspecification. To calculate 
Tnil), we also use Whittle estimator 'jf. under the model FAR(pk,d) for k = 1,2. 
We fix pi = p2 = 3 for n = 64 and pi = P2 = 5 for n = 128. 

We consider the following three alternatives, which admit the same forms as 
those in Hong (1996). Note that our alternatives are closer to the null hypothesis 
than those used in his paper. 

Alternative 1: Puvij) = 0.05 if j = and otherwise. 

Alternative 2: p^vU) = 0.05 if j = 0, sin(0.057rj)/(7ri) if 1 < J < 8 and 
otherwise. 

Alternative 3: Puvij) = 0.05 if j = 3 and otherwise. 

In the simulation of Alternative 3, we generate {ut, vt}t=-40oo,n by multiplying 
the square root of its covariance matrix, which is a sparse matrix with dimension 
(8002 + 2n) X (8002 + 2n), with a vector of (8002 + 2n) random numbers indepen- 
dently generated from standard normal distribution. We use the following three 
kernel functions: 

1. Bartlett (BAR) K{x) = 1 — if \x\ < 1 and otherwise; 

2. Tukey (TUK) K{x) = (1 + cos(7r2;))/2, if |x| < 1 and otherwise; 

3. Parzen (PAR) K{x) = 1 - Gx"^ + 6\x\^ , if\x\ < 1/2, 2{l-\x\f if 1/2 < |x| < 1 
and otherwise. 

For each kernel, we try three bandwidths: i?„ = [3n^'^], [3n'^'^] and [Sn*^'^], where 
[a] stands for the integer part of a. These rates lead to Bn = 6, 10, 15 for n = 64 
and Bn = 7, 12, 20 for n = 128. 

Tables [UU] show the proportion of 5000 replications in which the null hypothe- 
sis was rejected at 5% and 10% nominal significance levels. Both tests are carried 
out using upper tail critical values of the standard normal distribution. Note 
that the estimated standard error of the relative rejection frequency is given by 
^a(l -a)/5000, where a is the observed relative rejection frequency. Under 
model (jlip . it is seen from Table 1 that both tests are oversized at 5% and 10% 
levels with less size distortion for n = 128. The larger bandwidth and the Tukey 
kernel correspond to better size. Under model (|12p . the size distortion of Tn{6) is 
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noticeably larger than that of T„(7) uniformly in the kernel and bandwidth. We 
attribute this phenomenon to the model misspecification, which leads to inade- 
quate prewhitening for Tn{6). While for r„(7), the FAR(pfc,(i), A; = 1,2 models 
provide a decent approximation to the true data generating process, so the size is 
much less distorted. 

Tables [IH report size-adjusted power of our test statistics under three alter- 
natives. Under Alternative 1, the larger bandwidth corresponds to lower power. 
This is expected since the test statistics assigning weights to a large number of 
lags are less powerful than those assigning more weights to short lags in detecting 
the simultaneous cross-correlation of {ut} and {vt}. Under model (fTTIl . the test 
statistic Tn{9) seems slightly more powerful than T„(7) for all three alternatives. 
While for model (jl2p . T„,(7) outperforms Tn{&) in terms of power performance, 
which might be due to the model misspecification in forming Tn{6). Under Alter- 
native 3, the power for Bn = [3n^'^] is substantially lower than that for the other 
two bandwidths. Compared to the other two kernels, the Parzen kernel performs 
very poorly when Bn = [3n°'^]. Since the correlation of {ut} and {vt} occurs 
only at lag 3 under Alternative 3, any test that gives less weight to the 3rd lag 
squared cross-correlation tends to have less power. To illustrate this point, we 
note that when Bn = [3n'^'^], the (normalized) weights assigned to the 3rd lag for 
Bartlett, Tukey and Parzen kernels are 0.062, 0.056 and 0.019 for n = 64, 0.069, 
0.071 and 0.036 for n = 128. Since our results only require finite fourth moment 
of the innovation processes {ut} and {vt}, we repeat the above simulations with 
the innovations generated from a t distribution with 5 degrees of freedom. The 
performance of the tests (results not shown) is very close to that of the tests when 
the innovations are Gaussian. 

Since Tn{d) (r„(7)) is asymptotically equivalent to G„ (see ({2])), which is a 
positive linear combination of positive random variables, its distribution is right- 
skewed and non-normal in the small sample case. To correct for the size distortion 
observed in the simulation studies, we apply the power transformation method 
[Chen and Deo (2004b)] and compare the performance of the modified statistic 
with that of the original one (results not shown). It turns out that the power 
transformation method does not always lead to better size. For Hong's (1996) 
test statistics, similar findings have been reported in Chen and Deo (2004b), who 
attributed these to the fact that the mean and variance of the test statistics are 
based on the asymptotic argument, thus could be inaccurate in small samples. 
In our setting, a detailed check of the sampling distribution of our test statistics 
shows that this is indeed the case for certain combinations of the bandwidth and 
the kernel. To obtain more accurate approximation to the mean and variance 
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of our test statistics, one can follow the approach in Box and Pierce (1970) and 
Ling and Li (1997) to derive the asymptotic expression of nvar(p), where p = 
{Puv{—Bn), • • • , Puv{Bn)y ■ This is beyond the scope of this paper and will be an 
interesting topic for future research. 

In summary, we observe reasonable size and power properties for the proposed 
test statistics. The test statistics are oversized when the sample size n = 64 and 
128, and the size distortion becomes less severe when the sample size is larger, 
e.x. n = 400 (results not shown). Under three alternatives examined, the power is 
fairly high even for a moderate sample size n = 128 so long as the kernel and the 
bandwidth Bn are appropriately chosen. Further research on the selection of the 
bandwidth parameter Bn is needed and it seems to depend on the alternative under 
consideration. Also no particular kernel was found to outperform the other two in 
any situation considered. Empirically, model misspecification was found to yield 
inferior performance in terms of both size and power. A theoretical investigation 
of the impact of the model misspecification on the asymptotic distribution of Tn{6) 
is certainly of interest. We leave this topic for future work. 
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4 Appendix 

In the appendix, the constant C is generic and it may vary from line to line. The 
following decompositions ([13]) and (fT^ will be frequently used in the proof. Note 
that 

^ n / oo \ ^ oo 

«;x,(A) = --=j; J^a.^xi^, e**^ = A(A)u;„(A) + -^ 

™ t=i \j=o / y TTn ^.^Q 

where [/,>(A) = YJ^tZLj ^t^^'^ " ELi Since fx,x,{\) = (2^)-i|^(A)|2, we 

have 

...(A) _ MX)MX) ^^^^^^^^^r-oM^^\ (,3) 



V/xiX,(A) \AiX)\/V2^ " V^\A{\)\ 
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Similarly, we have 

..,(A) _ ^(A)^.(A) where i^.(A) = ^^-°y-i^^^"^ (14) 



and y,, „,(A) = X:r=w ^'te^*^ " EHi ^te''^- 

For the convenience of notation, write Aj = A{\j)^ Bj = B{Xj), Ii2j = 
IxiXii^j), Wuj = Wu{Xj), Wvj = Wv{Xj), Ruj = Rui^j) and Ryj = Rv{\j). Let 
gj = Aj\Aj\~^ and hj = Bj\Bj\~^. Further let Dk{X) = ELi^^**^- Denote by 
ay h = max(a, h) and a Ab = min(a, b). Define 

27r """^ 

fuvW = — ^ W{X - Xj)gjWuiXj)hjWy{Xj), 

" i=i 

27r:^' , 27r"-' 



/n«(A) = — X] - Xj)Iuu{Xj) and A™(A) = — ^ M^(A - Aj)I,,(Aj). 
Proof of Theorem l2.lt The conclusion follows from the following three assertions: 

n-l 

E (\fx,X,{Xl)\' - 47T^\fuv{Xl)\') = O^i^/K), 



1=0 



^ n—1 

-E/M.(Az)-l = Op(i?„-'/'), fc = l,2 



=0 

and z^;-o u^'^v i _n — jv O,l), 

which have been established in Lemmas 14.11 14.21 and 14.31 respectively. <C> 
Lemma 4.1. Under the assumptions of Theorem \2.1\ it follows that 

n-l 

E (l/xiX,(AOP - 47r2|A„(A0p) = Op{y^). (15) 
Proof of Lemma |4. II According to (|13p and (jl4p . we get 

fx^Xii^i) - 2vr/«,;(Az) = — ^ VF(Ai_j)[\/27r5ji(;t,ji?„j + V^mhjW^jRuj + RujRvj] 
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Then LHS (left hand side) of ([E]) = Gi„ + G2„ + G2n, where 



G 



In 



47r2 



«=0 

n-l 



4^2 _ 

— ^ ^ M^(Ai_j)[\/27rs'ji(;„ji?„j + \/2TrhjW^jRy,j + 
i=i 



To show Gin = Op{\/Bn), we note that Gin is smaller than 4(Gii„ + Gi2n + Gi3n) 
where 



Svr' 



Giin = ^ ^ W^(A^_J•)W^(Az_j/)5j^^;^^J•i^^j5^j/'^I;„J'-R, 



1=0 j,j'=i 



G 



12n 



— ^ ^ W{Xi^j)W{Xi^f)hjW—Rujhj,w^j'Ruj', 
" z=o i,i'=l 
n— 1 71—1 



47r2 



13n 



^ ^ W{Xi_j)W{Xi^y)RujRvjRuj'Rvj'- 



1=0 j,j'=l 



We shall only prove Gun = Op{-\/Bn), since the treatment for Gi2n and Gi^n is 
similar. Let A„(A^) = YJh=-B„ K'^{hbn)e^^^% where 6„ = 5,7^ In view of the fact 
that 



n-l 



Y,W{Xi.j)W{Xi.j,) = (47r2)-inA„(A 



(16) 



we have 



n-l 



Glln = 2'Kn ^ ^ gjgjiWujWuj'RvjRvj'^ni^j-j' 
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By Lemma 14.61 we get 



An 



2 n— 1 n— 1 



X] ^i'^uh'^uj[Wuj2Wuj!,MRvnRvj[Rvj2Rvj'J 

Sji 9j[ 9j2 9j'^ An ( Aj, _j J ) A„ ( ) 

^ n— 1 ^ n— 1 



j 1,32 = 1 



|c4(u) - 3| " ^ 



A^ ( „j,-j ) |E i?i,j2 



+j2-j[) 



ji,j'id2=i 

Let l(-) denotes the indicator function. The following fact 

|A„(A,)| = l(s = mod n)0{Bn) + 0{\Xs\~^)l{s / mod n), (17) 

which was stated in Equation A. 28 of Chen and Deo (2004), and Lemma [4.71 result 
in E|Gi„|2 = 0(log2 nBl/n^) + 0{Bl/n^/^) = o{Bn). Hence Gm = Op{^). 
To show G2n = Op{y/B„), we first show that 

47r^\/2 — n—i 

G21n ■■= ^ fuvi^l) ^ W{Xi^j)gjWujRvj = Opi^fSn). 

l=Q j=l 



n — 1 n— 1 



n 



Using p6j) again, we get 



87r^\/2^ 



21n 



J] J] T^(A;_,)1^(A,„ 



2ttV2it 



1=0 j,j'=l 



n 



■ ^ An{Xj'^j)gjgj'hjiWujWuj''Wyj'Ryj. 
j,j'=i 



So 



E|G: 



21n| 



n— 1 n— 1 



" ^ X] X] ^(^«Ji^«i(^«J2'^«i^)IE(^^t,jii?,;jil«t,j^-Ri;j2) 

X A„ ( Aj-/ )5fji ^/ij/ A„ ( Aj/ )gj;gj'^ hj,^ . 

Similar to E|Gii„p above, E|G2inP can be bounded by a sum of four terms in 
view of Lemma 14.61 For example, the first term corresponds to the case ji = j[, 
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J2 =j2, which is 

n-l [n/2] 

iij2=i ii,j2=i 

by Lemma 14.71 The bound for the other three terms can be estabhshed using 
Lemma 14.71 in a similar fashion. Hence G2in = Op{\/Bn). Following the same 
argument, the other two terms in G2n can be shown to be Op{\/B^). This completes 
the proof. 



Lemma 4.2. Under the assumptions of Theorem \2.1\ u'^YII^Zq fXkXi,{\) — 1 = 
Op{B-^'^), k = l,2. 

Proof of Lemma 14.21 We only deal with the case k = 1. Observing that 

n-l 

J2 W{Xi^j) = (27r)-^n for j = 1, 2, • • • , n - 1, 

1=0 

we derive 

71—1 „ n— 1 n~l r-ir/ \ \ T / \ \ -1 n—1 



1=0 1=0 j=i J^i^iy 31 



|2 



+V2TTgjWujRuj + V27rgjWujRuj] =■ Ji + J2 + h + Js- 

Let luuW = (2™)""'^! ~ ■iin)e**'*'P be the periodogram of the centered 

innovations ut — Un with Un = "-""^ Z^tLi ^t- Since luui^j) = luui^j) at j = 



1,2, ••• ,n - 1 and 27rn X]j=o -^""(^i) = z]"=i(^t ~ ""n) , we have J; 
n^"*^ Z]"=i(^t ~ "Un)^ = 1 + Op(n~"'^/^) by the central limit theorem. 
By Lemma |4.7|. we get 



27r 



var(J3) = — gjgj'COY{wujRuj,Wuj'Rv 



_ _ "-Ml' 

z — / 

— [cum(u;„j,ii„j,IZ}~7,i?„j,) + cov(7i;„j,II^ 

[n/2] 

+cov{wuj,Ruj')coY{w~7,R^j)] = 0{n~'^) ^ {jj'r^^'^ = 0(n~^). 
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A similar argument yields var(J2) = 0{n ^log^n). The conclusion follows. 



Lemma 4.3. Under the assumptions of Theorem \2.1{ we have that 



^/2Bnd{K) 



iV(0,l). 



Proof of Lemma [43) Let Gjji := gjhjgj'hji and Hn := 47r^ SILq^ \fuvi^l)\'^- Then 
by dUD, 



n 

1=0 j,j'=l 



An 



^ Gjj>An{Xj-j')WujWyjWuj'Wyj' (18) 



n 



^ n n n n n—1 
tl=lt2 = lt3 = li4 = l ij'=l 

Denote by 5n(A:i,fc2) = EiJ'Li G,j'A„(Aj„jv)e*fci^^e-''=2V and <7„(A;) = gn{k,k). 
Note that gn{k) is real. Write 

^ n-l (n+fci)An (n+fc2)An 

-f^n = — ^ ^ ^ Ut^Ut^Vt^^kiVt2~k29n{h,k2) 

fci,A:2=l-n ti={l+fci)Vl t2={l+fc2)Vl 
^ n_l (n+fc)An {n+fc)An 
= -f^On + — EE E Ut-^Ut2Vti-kVt2-kgn{k) = Hon + Hin + H2n, 



n 

k=l~n ti = (l+fc)Vl i2 = (l+fc)Vl 



where 

Hon 



^ n—l (n+fci)An (n+fc2)An 

-3 E E E UtiUt2Vt^-kjVt2-k29nikl,k2), 



fcl^fc2=l-nii = (l+fei)Vl t2={l+fc2)Vl 
71—1 n 



Hln = -3E E Ut^Ut2Vt^-kVt2-k9n{k), 
k=0 ti,t2=l+k 
^ 71—1 n 

H2n = — E E '^hVt2Ut^_kUt2-k9n{-k). 



k=l ti,t2=k+l 
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We shall first prove -ffon = Op{-\/Bn). Note that 

n-l n-1 (n+A;i)An (n+fc2)An (n+fei)An (n+fc^)An 

var(i^o„) = H 

ki^k2=l-n k{j^kL,=l-n ti=(l+fci)vl t2 = (l+fc2)Vl t{ = {l+k{)\/l i^={l+^)Vl 

E [uti u^/^ ]E [ut J vj/ v^/^ __fe^]fi(„ {ki,k2)gn{k'i,k'^). 

By Lemma 14.61 there are at least four restrictions on ti , t2 , t'l , ^2 iki,k2, k[ , for 
each non- vanishing term in var(ffon)- For example, one of such terms, denoted by 
Jin, corresponds to ti = t2,t[ = t'2, ki = k[, k2 = k'2, i.e. 

n-l (n+ki)A{n+k2)An 

Jm = ©(n-*^) Yl Y \9n{ki,k2)\^ 

ki^k2=l-nti,t{={l+ki)V{l+k2)Vl 
n-l 

= 0{n-') Yl \9n{ki,k2)\' 

ki ,k2=l—n 
n-l 

= 0{n-') Y \^ni\n^,0An{Xj2~r2)\ 

In view of (HZ]) and the fact that Efc=i-n = l(s = mod n){2n - 1) - l(s / 
mod n) for s G Z, we have Ji„ = 0(log^ n) = o{Bn) under Assumption l2.31 Other 
terms in var(i/on) can be bounded by O(log^n) in a similar way, so var(i7on) = 
o{Bn) and H^n = Op{y/B^). 

Write Hin = Hun + -H'i2n and i?2n = -f^2in + H22n, where 

2 n— 1 n ti — 1 

-fflln = "sX] X] X] '^ti^t2'"tl-kVt2-kgn{k), 
k=0ti=k+2t2=k+l 
2 n— 1 n ti — 1 
H21n = ~Y Y Y '"tiVt2Ut-^-kUt2-kgn{-k), 
k=l ti=k+2t2=k+l 
^ n—1 n ^ n-l n 

Hl2n = ■:^YY ^tVt-k9nik), H22n = :^Y Y ^tUt-k9ni-k) ■ 
k=Ot=k+l k=lt=k+l 

We shall show that var(//i2n) = o{Bn) and \ai{H22n) = o{Bn)- Note that 

^ n— 1 n n 

var(iJ22n) = -5 Y Y Y 9ni-ki)gni-k2)cov{vf^uf^_k^,vluf^_k^). 

fcl,fc2 = l tl=fcl + l t2=k2 + l 



n-l 



fci,fc2 = l- 
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Since cov(7;4n4_;^^,<n4_fcJ = E(t>,>4)E( 



1 = [c^iv) + 2]l{ti 



t2,ki / A:2) + [c4(M) + 2]l(ti-fci =t2-k2,ti /t2) + [(c4(?^)+3)(c4(M) + 3)-l]l(ti = 

t2,ki = k2), we can write var(i?22n) into a sum of three terms, which correspond 
to ti = t2, ti — ki = t2 — k2 and ti = t2,ki = k2- When ti = t2, the corresponding 
term is 



J2n = 0(n-6) Yl E 5n(-^-lK(-fc2) = 0(n-6) ^ 

fcl,fc2 = l il = (fci+l)V(fe2+l) il J2,J2 = 1 

n— 1 n 



|An,(A._,;)||A„(A,,_,.)| 



n-1 



0{Bln-') 

n-2 



fci,fe2=l ti = (fci+l)V(A:2+l) 

^ Ai-i(Aj,_jj)Ai-i(Aj2_j 



J2- 



mi ,1712=2— n 



ii=2 
n 



^ Ai-i(A„,jAi-i(A 

ti=2 



m2. 



Since Dk{Xs) = 0(|A,|^i)l(s / 0) + fcl(s = 0) for s = 2 - n, 3 - n, • • • , re- 3, n-2, 
under Assumption 12.31 we have 

n-2 n-2 

|J2n| = 0{Bl/n') Y (An^,A^2)-'+0(5nV) E ^-1 
mi,m2=l m.i=l 
n 

+0(B2/n4) ^(ti _ 1)2 = 0{Bl log^n/n) = o{Bn). 

ii=2 

Similarly, we can show that the other two terms corresponding to ti — ki = t2 — k2 
and ti = t2, ki = k2 are of order o{Bn)- So -ff22n = Op{-\/Bn)- Regarding -ffi2n; we 
note that i?i2„ = re'^ YJi=i Ut'^t9n{0) + Ylk=i StLfc+i '^t'^t-k9n{k), where the 
variance of the latter term is of order o{Bn) by the same argument as above and 
the former term is easily shown to be Op{y/B^). Therefore Hn = Hun + i?2in + 

Op{^/B^). 

Since E(i?„) = Bns{K) + 0(1) and var(i7„) = 2Bnd{K){l + o(l)), as proved 
in Lemma |4.5[ our conclusion holds if we can show 



Hlln + H2ln 



A^(0, 1), where al = 25„d(K). 



(19) 
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INOte inat l^k=l l^t=k+2 l^s=k+l — l^t=3 2^s=2 l^k=l ^^"^ l^k=0 2^t=k+2 l^s=k+l — 

Et=2 El=\ Efc=o- We can write Hun = Et=2 Wmt and H^m = Zts ^2nf , 
where 



2 t—i s—i 2 *~i 



= — ^^Mtn^ft_fc?;s_fc5'„(A;) and W2nt = —'^^vtVsUt-kUs-k9n{-k). 

s=l k=0 s=2 k=l 

Then Wint and VF2nf form martingale differences with respect to Tt, where J^t is the 
cj-field generated by {us,Vs}1=_oq- Letting Wnt = Wint+W2nt, then Hiin+H2in = 

^^^^ _|_ n~^Win2i where the latter term is easily seen to be Op(l). 

We shall apply the martingale central limit theorem of Brown (1971). It suffices 
to verify the following two conditions: 

n 

a~''n~^Y,Wl^pl, where W^^ = nwl\Tt-i) (20) 

i=3 
n 

and a-'^n-'^^W.[W'^^l{\Wnt\ > eracj„}] 0, for any e > 0. (21) 

Since cr-'^n^'^ E"=3lE(W^4) = cr-2var(iJii„, + if2i„) + o(l) = 1 + o(l), as we have 
shown, (j20p is implied by 

var ^a-2n-2 ^ t|/2^^^ = ^(1), where Wf^, = E{W^JJ^t-i), 3 = 1, 2, (22) 

which is established in Lemma |4.4[ To prove (|2ip . it suffices in view of (a + b)^ < 
8{a^ + b^),a,b£R to verify 

n 

a-^n-^Y.^{Wf^,) = o{l), i = l,2, 

which can be shown by a similar and slightly simpler argument as in the proof of 
Lemma |4.4[ The details are omitted. 

Thus ([19]) is established and the conclusion follows. 

Lemma 4.4. Under the assumptions of Theorem \2.1l the random variable Wj^^ 
defined in satisfies 



var c7;^2^"2 ^ Wf^t = 0(1), j = 1, 2. 



t=3 
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Proof of Lemma [4.4l We shall only show the case j = 1, as the treatment for j = 2 
is similar. Note that 



't-l s-l 



t-l 



Xl UtUsVt-kVs-kgn{k) + UtVt ^ UsVsgn{0) \^t-l 
\s=2 k=l s=l J 

't-l s-1 \ 2 /t-l \ 

^"^UsVt-kVs-kgnik) + ^UsVsgn{0) 



vs=2 k=l / 

Let Jo := var(X;"=3(Es=i^«"^ffn(0))^) and Ji := var(X;"=3(Es=2 Efc=l ^s^^t-fc 
V s-kgnik))^) ■ The conclusion follows from = o(cr^7i^^), A; = 0, 1. We consider 
Jl first. 

n I tl — 1 si — 1 S2 — 1 

Jl = X] X] X] X] Us-,Us2Vt^-kiVs-,-kiVt-i^-k2Vs2-k2gn{kl)gn{k2), 

tl,ts=3 \ Si,S2=2 fei = l ^2 = 1 

ts-l Ss-1 S4-1 

X] X] X] '^SiUs^Vt^-kiVs3-k3Vt3^kiVsi-k4gn{k3)gn{k4) 
S3,S4=2 fc3 = l /C4 = l 

?1 tl — 1 Si— 1 S2 — 1 ta — 1 S3 — 1 S4 — 1 

= X] XIX^X^X^X^X] 9n{ki)gn{k2)gn{k2,)gn{ki)Jii, where 

tl,t3=3 si,S2=2 fei=l k2=l S3,S4=2 fc3=l A;4=l 

J^l = CO\{Us^Us2Vt^-k^Vs-^_k^Vt^-k2Vs2-k2^'Us.^Us4Vt.^-k3Vs3-k-iVt3-kiVsi-ki) 
= E{Us-,Us2Us._,Us^)K{Vt^_k^Vs-^_kiVt^-k2Vs2-k2Vt3-k3Vs3-k3Vt3-k4Vs4-ki) 
— fci ^Si — fci ''^tl — fc2 ^S2 — fc2 )E('Us3 ^S4 ) 
XE{Vt.^_k3Vs3-k3Vt3-k4Vs4-k4) = -^111 -^112 " Jll3- 

In the above expression, Jm = E.{usj^Us2Us3Us_^) = l(si = S2,ss = S4) + l(si = 

S3, S2 = S4) + = S4, S2 = S3) + {Ci{u) - 3)l(si = S2 = S3 = S4), Jll3 = = 

S2, S3 = Si, ki = k2, k^ = ki) and 

J112 = nvt 1— fci'i'si—fci 'f^ti—fc2^S2—A:2 ^43-^3 ^S3—A;3 '^t3—k4 ^34 — ^4) 
= ^ cm'n(i;i^. , € 51) • • • cum(i;i^. , ij G 5^), 

where is over all partitions g = {gi U • • • U gp} of the index set {ti — ki,ti — 
^2,^3 — ^3; ^3 — ki, si — ki, S2 — k2, S3 — k^, S4 — ki}. Since E(i)j) = 0, only partitions 
g with ^gi > 1 for all i contribute. We shall divide all contributing partitions into 
the following several types and treat them one by one. 
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1- 7^51 = #52 = 4. A typical term that contributes is 

cum{vt^_ki , Vtj^_k2^'"t3-k3^Vtg~k4)cnm{vs^-ki,Vs2-k2^Vs3-k3,Vs4-k4), 

which equals to l{ki = k2, = k^, si = S2, S3 = S4, ti — ki = t^ — k^, Si — ki = 
S3 - kz). 

2- #51 = #52 = 3, #53 = 2. A typical term is 

cum(Dj^_fci , fti_fc2, Ut3_fc3)cum(-Us-^_fc^ , Vs2-k2,vt3~k4,)co-v{vs^_^-k3,Vsi^k4), 
which is l(/ci = k2-, si = S2,ti—ki = t^—k^, ^3— ^4 = si—ki, s^—k^ = 54—^4). 

3- #51 = #52 = #53 = #54 = 2. One such term is 

COvivt^_ki,Vt^-k2)cOv{vt^-k3,Vt3-k4)cOv{Vs-^-k^,Vs2~k2)cOv{Vs3-k3,Vs4-k4),{'^^) 

which is l{ki = k2, ^3 = k^, si = S2, S3 = S4). Note that (f23]) multiplied with 
l(si = S2, S3 = S4) in Jill cancels out J113. Thus all non- vanishing terms in 
Jii involve at least five restrictions on the indices ti, ^3, si, S2, S3, S4, ki, k2, k^, k^. 

In the following we will find a bound for only one such term since the derivation 
for other terms is similar. For example, one of the terms in Ji, denoted by E* , 
corresponds to the case /ci = A;2, A:3 = A:4, si = S2 = S3 = S4, i.e. 

n (ti-l)A(i3-l) si-1 

^* = E E E 9l{ki)9l{k3). 

ti,t3=3 si=2 ki,k3=l 

Recall gn{k) = TJj7jLiGjj'K{\j~j')e'''^^-^'. We have 



n— 1 n— 1 

n. .,n. ., A (\ . (\ - , , 

^]2-]2> 



E* - ^ Gj3j/Gj4j/A„(Aj^_j/)A„(A 

h ,32 ,j2 = l i3 J3 ,j4,i4=l 

n (ii-l)A(t3-l) 

A„(A,3_,,)A„(A,,_,,) E E E e^'^'--l+--^e^^^^--3+^4-4. 

ti,t3=3 si=2 fci,fe3=l 

Then 

n-2 

\E*\<Cn'' Yl |A„(A^J||A„(A„2)||A„(A^3)||A„(A^J| (24) 

mi ,m2 ,m3 ,m4=2—n 
n {ti-l)A(i3-l) si-1 

^11^3=3 Si=2 ki^k^ = l 
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_. TP* \ 

• ^■mim2m3m4/' 



It is not hard to see that E^^^^^,^^^ is 0{\Xmi+m2\ ^lAmg+m*! '^n^) if mi + m2 ^ 
{0, ±n} and 7713 + 7714 {0, ztn}; 0(|Am3+m4l~"'^?^^) if mi + m2 G {0, ±n} and 
"^3+772.4 ^ {0, ±?i}; 0(|Ami+m2 1""'^'^^) if mi+m2 ^ {0, ±71} and m^+nii G {0, ±?i}; 
0(71^) if TTii + 7712 G {0, ±?i} and m3 + m4 G {0, ±71}. 

Combined with pT|) . we can derive from ([Ml) that = 0(77^^S^ log^ 77) = 
o{n^'^B^). Therefore | Ji| = o(cj^77-^^) and a similar argument yields Jq = o{a^n^'^). 
The conclusion follows. 



Lemma 4.5. Under the assumptions of Theorem \2.1\ the random variable Hn in 
/ flgj) satisfies 

E{Hn) = Bns{K) + 0(1) OTid var(F„) = 2Bnd{K) + o(S„). 

Proof of Lemma l4.5t Recall from Lemma 14.31 that Hn = Air'^n"^ j'Li ^n(Aj_j') 
Gjj'WujW^-w^w^j'. Since E(7x;uj^I^) = E(mi^Wyj') = {2iT)~^l{j—j' = OmodTi), 
we have E{Hn) = (1 - 1/n) T.h^=~Bn ^^(hiK) = Bns{K) + 0(1) under Assump- 
tion O 

We proceed to calculate var(if„,). Note that 
var(if„) = -|- J]] Gj^jjG~~7A„(Aj-^_j/)A„(A^2_j-^)C(ji,i2,il,j2), 

where by Lemma 14.61 

C (ji , j2 , i'l , J2 ) := cov(uj„j^ u)— • W—7w^,j'^ , uh^Wyj^ ■ Wuf^w—r) (25) 
= E(7/;t,j^ u)„j/ u;„j2 -w^jv )E{w~w^ji_^ ■Wyj^W~r) 
-(167r4)-4(ii = j0l(i2 = j'2) = (167r4)-ni(ji +i2 = n)l{j[ + j2 = n) 
+l(ii = J2)1(J2 =j'i) + (C4(n) - 3)(c4(i;) - 3)n~H{ji + - j[ - j2 = 0, ±n) 
+21{h+j'2 = n)l{h = j[)l{j2 = 32) + 21(ji = r7)l(ji = i2)l(j2 = Ji) 
+21(ji = j2 = j[ = j'2) + (C4(n) + c^{v) - 6)77"Hl(ji = Ji)l(j2 = ia) 

+ l(il = 77)l(i; +i2 = n) + = J2)l(i2 = fl)}]- 

A simple calculation shows that only the first two terms in C{ji, j2, j'l, j'2), i-e. 
(167r4)-i[l(ji + = r7)l(j; +i2 = n) + l(ii = j2)l(i^ = jO] contribute to the 
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dominant term. Therefore, 



var (Hn) 



2(l + o(l)) 



71-1 

n-l Bn 



jl,jj = l hl,h3 = -Bn 

^2""^^^^ E K'ihibn)K\hsbn) E ^ cos(A:A;,,„,,3) 

hi,hs=—B„ k=l 



4(l + o(l)) ^ , , 



fc=i 



where an(/c) = Eh=-B ^^(^^n) cos(feAft) and = Ylh=-B K^i^bn) sm{kX 



h)- 



By Lemma 4 in Wu and Shao (2007), 



k=l 



Bn 



.h=l 



+ 0{nBl). 



Since K{-) is symmetric, bn{k) = 0. Thus 



h.=l 

This completes the proof. 



Proof of Theorem 12. 2t The proof fohows the argument in the proof of Chen and 
Deo's (2004) Theorem 4. Let ^0 = (^io'^20)' and 6 = (^'1,^2)'- It suffices to show 
that n{Tn — T„,(^)) = Op{\/Bn), which follows from the following two assertions: 



n-l 



Y.(\fxiX,{Xi)\'-\fx,x,{Xi)\') = op(v^), 



/=0 

2tt 



n-l 



n 



E{/m(Az)-/m(Az)} = Op{B-^/^), k = l,2. 



(26) 

(27) 



1=0 



We shall only provide a proof for (j26p . since the treatment of (j27p is similar. 
For e = {e'M let G{\j,\j,-e) = Ul^J-^^\Xf,0k)f;:^^\Xjr,dk). Further let 
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1*12 j = ^i2j/V/i(Aj;^io)/2(Aj;^2o). Then the LHS of m is 
47r 



E W^(A«-i)VF(V,v)/i2,/i2/{G'(A,-,A,v;^)-G(A„A,v;^o)} 

«=0 j,j'=l 



_ 1 A /X ■7;^ {G'(Aj, Aj/;g) - G(Aj, Aj/;go)} 

- - 2^ ^ni^J-/)A2,A2i' G(A,,A,.;^o) ' 

where we have apphed (|16p. Let 7i and /2 be the index set of Qi C M"^^ and 
02 C M''^ respectively. For every \j and Aj', a Taylor series expansion yields 

G{Xj,\jr,e) - G{Xj,Xjr,eo) = y^i^iu -^lOu) — ^ ia ^ ' 



:\^ (n n ^ OG{Xj,Xj';tfo) \ - d G[Xj,Ajr,(fjj,) ^ 



where ^jj/ = {dijji,62jjiy = + otjj'i^ ~ ^o) for some ajji E [0, 1]. 

ForA; = 1,2, let gku{^;Ok) = d\og fk{\]6k) / dOk^ and gkuv{>^]Ok) = d"^ / dek^dOk^ 
log/fc(A;^fc), u,v £ Ik- Then for u G 4, 

^ q1 ' = --^G{Xj, Xj';9o)[gku{Xj;6ko) + 9ku{^j'',9ko)]- (28) 

o^k^ 2 

Let Auv{Xj, Xj']djj') be the (ii,t')th element of the matrix d'^G{Xj, Xj']djj')/dd'^ . 
Then 

^™(Aj, Aj/;^j,,v) = -G{Xj, Xjr;djj>){[gku{Xj;6kjj') + gku{^j'',0kjj')]'^ 

-^gkuvi^j'jSkjj') + 9kuv{^j'',0kjj')]}, U,V e Ik,k = 1,2. 

= ^^(-^i' A/;^jj'){bfci«(Aj;^feijj') + 5fci«(Aj';^fcijj')][9fc2i>(Aj;^fc2ji') 
+5fc2t^(Aj';^fc2jy)]}> ''^ e e 42, (^1,^2) = (1,2) or (2,1). 

Since \0k — dko\ = Op{n~^/^) and (j28p . (j26p follows from the following two state- 
ments: 

1. For each n G 4, A; = 1, 2, 

L„ := E An(Aj_jv)/*2j7^[5fc„(Aj;^fco) +5fcn(Ajv;^fco)] = Op{n^/^ y^). (29) 
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2. For every (u, v), 

Mn := J2 An(A,-,0/i%-7^ ^g^^^'^f.g = (30) 



■ ., 1 G{Xj,Xj']6o) 

By (ini) and L n can be written as a sum of 16 terms, of which the dominant 
term is given by 

n-1 

Lin := 47r^ ^ An{Xj-j')Gjj'WujW;jj-w^Wjjj>[gkii{Xj]0ko)+gkuiXj/-,0ko)]- 
Then E|Li„|2 is 

n— 1 n—l 

IGtt^ ^ ^ An{Xj^_j>^)An{Xj^_j0Gj^j'^Gj2j>^[gku{Xj^]eko) + gku{>^j[',Oko)] 
ii>i(=ii2j2=i 

xbfc«(Aj2;5fco) + 5'fcn(Aj^;^fco)]IE(u;„jiIZ}— 7 • w;;j;w^j;^)E{w;—w^j>^w^j2W-r). 

In view of Lemma 14.61 and (I25p . the dominant term of E|Li„p is a sum of three 
terms, which correspond to ji = j[ and j2 = Ji + J2 ~ ''^ ii + J2 = ''^ as 
well as ji = j2 and = i^. Since \gku{^]0ko)\ < C|A|~'' under Assumption 12.71 
all the three terms can be bounded by 

n-1 

CBl i^n^nr'' = 0{n^Bl), for 5 G (0, 1/2). 

jlj2 = l 

For the other terms involved in E|Li„|^, a tighter bound than 0{'n?Bn) can be 
established using Lemma O and (HZD. So Li„ = Op{nBn) = Op(n^l'^B]l'^). The 
remaining 15 terms in L„ can be shown to be Op{n^^'^B}/'^) using the bounds estab- 
lished in Lemma 14.71 Since the proof does not involve additional methodological 
difficulties, we omit the details. 

Now we prove (pOj) . Note that Okjj' = iPkjj'^dkjj')' , k = 1,2. Under Assump- 
tions [2]6] and [221 we have 

\AUXj,X,r,ej,>)\ < C{X,Xj,f^n'+^2n'^X-^ + X-,^f 
uniformly in (u, v). By Lemma 14.71 and ()13p . uniformly in j = 1, 2, • • • , n — 1, 

^l^a^'^i^^^ ^ 2[2.E|^.„,f +E|i?„,f ] < a 
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Hence, by the Cauchy-Schwarz inequality, uniformly in j, j' = 1,2, ■ ■ ■ ,n — 1, 

E{/x,x,(A,)} E{Ix,x,iXj)} 



m2jihf\< sup f(l p \ sup fiTft \ <c- (31) 

In the sequel, we shall treat the following several cases separately. 

1. di > diQ and ^2 > ^20- Then dijji > diQ and d2jj' > ^20 for all which 
implies that \G~^{Xj, Xjr,6Q)Auv{Xj, Xj';6jji)\ < C{X~^ + Xj,^)'^ . In this case, 
we get 

n-l 

|M„|l(di > dl0,d2 > d20) <C lA*2i^l|An(A,_,v)|(AT^ + a/)2. 

Then by and ([I7D, E|M„|l((ii > (iio,(i2 > (^20) is bounded by 

n-l 

C 5^ |A„(A,_,OI(A-' + = OinHogn). 

So |M„|l((ii > diQ,d2 > d2o) = Op{n?\JBn) under Assumption 12.31 

2. di < dio and(i2 < c^20- Denote by A^.-^ (Xj) = fl^'^{Xj;6ko)fi^^^'^iXj;dkjj')-l, 
k = 1,2. Under Assumption 12.61 and ([7]), for j = 1, 2, • • • , n — 1, 

|Af (A,)| < C\9k-eko\X-^'''"^^\ k = l,2. (32) 

Let Njji := G ^{Xj, Xji',d jji)A^y(yXj, Xji',0 jji^ . Then 

n-l 

\Mn\l{di < dio,d2 < d20) < \^n{Xj-j')\\Iujl!^\\Njj'\ 

x(l + |Af (A,-)|)(l + |Af (A,0|)(1 + |Af (A,)|)(l + |Af (A,-0|). (33) 
The RHS (right hand side) of (|33|) consists of 16 terms, one of which is 

n-l 

J„:= Y \K{X,.f)\\It2j7!^\\^{'\Xj)m,y\. 
According to §2) and the fact that \Njj>\ < C{X~^ + X~,^f, we have 

{n-l 
Y |An(A,-_,-OI|/r2,7S;:|A-(''-+'^^(A-^ + A/)2 
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where the expectation of the bracketed term above is bounded by 

n-l 

C |An(A,_,')|AT('^-+'^)(AT^ + A/)2 = 0(n2logn) 

for G (0, (l-dio-K)/2). Together with the fact that |^i-^io| = Op(n-^/2), 
we get Jn = Opin^yJBn)- The other terms on the RHS of (j33p can be shown 
to be Op{n'^\/Bn) in a similar fashion. Hence |Mn|l(di < (^10,^2 < c^2o) = 

3. di > dio,d2 < d2o and d2 > ^20, di < c^io- These two cases can be handled in 
a similar manner as in the previous two cases and the details are omitted. 

The proof is now complete. 



Proof of Theorem 12. 3t Similar to the proof of Theorem 12.21 it suffices to show the 
following two assertions: 

n-l 

Y.(\fx^x,{Xi)\'-\fx,xM\') = op(v^), (34) 
1=0 

-E{/m(AO-/x,x,(AO} = Op{B-'/% k = 1,2. (35) 

1=0 



We only prove (j34j) . as the proof of (j35|) is similar. Let 

=nLi/x;£(A,)/x;£(A,'), H,,,{f)=uLj-'/\x,;j,)f-'/\x,r,^,). 
Denote by 70 = (7'io>72o)'- Applying ([16|), the LHS of ([MI) is 

1=0 j,j'=i 



^ n— 1 

- E An(Aj_j')/l2j/l2j'{[^fij'(7o) ~ ^jj'] + [-f^ji'(7) ~ ^jj'ho)]} = 'S'ln + 'S'2n- 



n 
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Note that Assumption 1 2 . 8 1 implies that \Hjji{'yQ)H ~ 1| C'X]fc=i Ylu^pf,+i 
Therefore, bv (fT7|) . (I3TD and (l9|). 

n-l 

nSin\ < n-' \Ar.{X,^,,)\m,^7l^\\Hj,,{j,)H-^-l\ 

2 oo 
k=l /=pj;+1 

Regarding 52^, we have 

The remaining proof largely follows the argument in the proof of Theorem 12.21 
Here we only outline key steps. Let = ■ ■ ^pk + 1} be the index set for 7^^.. 
For each we apply a Taylor's expansion and obtain 

H,Al) - H.A-fo) = EE(^^^-^^oJ^^f^ (36) 

fc=l liG Jfc " 

where jjj, = {I'ljj' ,l'2jj>)' = 7o + "ii'(7 -7o) for some a^j/ G [0, 1]. 

Let hkuiX;Jk) = <91og/fc(A;7fc)/57fc„, hkuv{klk) = log /fc(A;7fc)/57fc^,^7fc. 
for u,v £ Jk and Buvi^j, be the (u, t;)th element of the matrix d'^Hjji{')jj,)/d')'^. 

Then we have 



dHjj'jjo) ^ _1 
dlk^ 2 



^^jj'(7o)[^fc«(Aj;7fco) + hku{Xj';i/ko)] 



and 



Buv{>^jAj']ljj') = -^Hjj>{'yjj,){[hku{Xj;Jkjj') + hkui>^j';i/kjj')f 

-2[hkuv{>^j;ikjj') + hkuv{>^j';ikjj')]}^ u,v e Jk,k = 1,2. 

= -^Hjj'i'yjj'){[f^kiu{>^j;'Jkijj') + ^^n«(Aj';7fcij/)][/ifc2t'(Ai;7fc2ii') 

+^fc2f(V;7fc2jy)]}>^ G Jfci,^ G Jk2,{ki,k2) = (1,2) or (2,1). 

Note that Hjji{')Q)H~^ < C under Assumption 12. 8i To prove S2n = Op[-\/B^), it 
suffices in view of Assumption 12.91 and p6p to show that 
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1. Uniformly in n G J^, k = 1,2, 

2. Uniformly in u G Jk^ and v G Jfcj, A;i, A;2 = 1, 2, 

j jj jj \ '0/ 

It is easy to see that Assumption 12.81 implies that uniformly in A, {u,v) and 7^, 

\hku{X;ik)\ < C\X\~\ |/ife™(A;7fc)| < C\X\-^ for any 6 > 0. 

Similar to the treatment of L„ in the proof of Theorem [221 we can derive E|L* p = 
0{n'^B^) = o{'n?Bn/p\) under ([10]). A probabilistic bound for M* can be estab- 
lished in the same way as that for M„, (see (|3Up ). Here we omit the details but 
mention the following fact (|37|) . which is needed in the proof. Under Assump- 
tion [2]8l we have 

l/r^(A;7<") - /,!'^(A;7f )l < CI7L" -7f l/r(A;7f ) (37) 

uniformly for all A and all 7^,^^ and 7^,^^ such that d^^^ < d^^^ . 

The conclusion is established. 

Lemma 4.6. For any ti,t2,t^, G Z, E[utj^ut2Ut3Ut4] = (c4(m) — 3)l(ti =^2=^3 = 
U) + l{ti = t2)l{t3 = U) + l{ti = t3)l{t2 = U) + l{ti = U)l{t2 = h). Further, 
for 1 < ji,j'i,j2,j2 < n-1, E{wuj^w—7Wuj2W~l) = (47r2)-i[l(ji = j[)l{j2 = 
32) + + 32 = n)l{j[ + = n) + = i^)l(i2 = j'l) + (C4(u) - 3)n-il(ji + 

i2-ii-i2 = o,±n)]. 

Proof of Lemma 14.61 Since {ut} are mean zero iid random variables, the first 
assertion follows from the following fact: 

¥.[ut^ut,,ut.^ut^] = cum{ut^,ut2,ut.^,ut^) + co\{ut^,ut.^)co\{ut.j,,Uti) 

+cov(Mti 1 ut.^)co\{ut2 ,uti)+ cov {ut^ , ut4)cov(ut2 ,Uts). 
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Applying the first assertion, we have 

1 



n —t2X., +tsXj„—t4X./ ) 



<l,i2,t3,t4 = l 



1 



n 



ti,i2=l 
n 



-it2X-i , ;/ 



il,i2 = l tl=l 



Since Y17=i ^ ~ '^^(j — ™od n), the second assertion follows. 



We now provide an auxiliary lemma on the bound of the second and fourth 
cumulants of Ruj ■ 

Lemma 4.7. Under Assumption\2.4\ we have that 



\cOv{Wuj^,Ruj2)\ 


= 0{\32 


-1/2) 




\cov {Ruji, Ruj2)\ 


= o{\h 




-1/2) 


\^{Ruji Ruj2^ujs Ruj4)\ 


= 0{\h 






c\xni(^Wuji J Ruj2 ) '^ujs ) -^Mj^) 1 


= o{\n 




-1/2) 


1 "^{Wuji Ruj2 W^uj-i Ruj 4 ) 1 


= 0{\j2 




-1/2) 



hold uniformly mji,j2,i3,j4 = -[n/2],--- , -2, -1, 1, 2, • • • ,[n/2]. 

Proof of Lemma 1121 Let Kn{\) = |i:>„(A)| V(2vrn) be Fejer's kernel. We shah first 
state a useful fact: 

■= r ^TT - 1 ^"(^ - ^i)^^ = O(Ur'), as n ^ oo (38) 

holds uniformly in j = —[n/2], • • • , —2, —1, 1, 2, • • • , [n/2]. The proof of ([38]) is ba- 
sically a repetition of the argument in Robinson's (1995) Lemma 3 and is skipped. 
Note that 

^iRujiRuj2^uj^Ruj4) — cum(^ujj , Ruj^, Ruj-^j Ruj4) ~\~ cov(/2^j-^, R^^j^) x 
cav{Ruj.^,Ruj^) + cav{Ruj^, Ruj.^)cav{Ruj2, Ruj4) + cov{Ruji, Ruj4)cov{Ruj2, Rujs) 
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and 



COv{Wuj3 , Ruji) + COv{Wuj^ , Wuj.^)cO\{Ruj2 , Ruji) + COv{Wuj^ , Ruj4,)cOv{Ruj2, Wnj-^) ■ 

Let Tj{X) = [A{—X) — Aj]\Aj\~^ . We shall find a bound for each term on the RHS 
of the equations above. After some straightforward calculations, we have 

1 



\C0V {Wuj-,,Ruj2] 



2'Kn 



\cov{Ruji,R 



■UJ2J 



< V27r 

1 

n 



D{Xj, - X)D{\ + Xj^)Tj,{X)dX 

1/2 



Kn{X,, - X)dX 



n 



1/2 



0(i2 



1/2^ 



and 



Tj,{X)T,,{-X)Dn{X + Xj,)Dn{-X + X,,)dX 



where we have applied the Cauchy-Schwarz inequality. Let = [— 7r,7r]^. By the 
Cauchy-Schwarz inequality and the periodicity, 

|C4('U)| 



I cu.in.(^Rujj , Ruj2 ) Rujg ) Ruj4 J 



[ Dn{X,, - Xi - X2 - X3)Dn{Xi + X,,) 



Dn{X2 + Xj,)Dn{X^ + A,jT,-,(-(Ai + A2 + X^))Tj,{Xi)Tj,{X2)Tj,{X^)dXidX2dX^\ 



[ Kn{X,, -X1-X2- X^)\Tj,{-{Xi + X2 + X3))\^dXidX2dX3 

/ Kn{X, + A,,)K„(A2 + X,,)Kn{X3 + A,' J jT,, (Ai ) ^ |T,3 ( A2) ^ X 



1/2 



\T,,ix-s)\'dx,dX2dx,Y/' < cp];^p]!^pli'p]'^ = o (ntij, 



-1/2 



Similarly, 



I CWUliWiiji , Ruj2 1 ^ujs 5 Ru]i ) 



ca{u) 



27rn2 



/ Dn{-{Xi + X2 + X:i))Dn{Xj, + Xi] 



Dn{X2)Dn{Xj, + A3)T,,(Ai)T,,(A3)(iAi(iA2(iA3| 



< C 



n3 



KniXj, + Ai)K„(Ai)K„(A,, + A3)|r,,(Ai)|2|r,,(A3)pa!AidA2dA3 



1/2 



j^^ Kn{Xi + X2 + Xs)dXidX2dX3 < CP]1''P]1^ = 0{j-^'''jl^'^' 



Therefore the conclusion follows. 
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Table 1: Rejection rates in percentage under the null hypothesis: (a) when the data 
are generated from model pT]) . (b) when the data are generated from model (fT2|) . 



(a) 








Tnie) 




Tni'l) 


n 


Bn 


a% 


BAR 


TUK 


PAR 


BAR 


TUK 


PAR 


64 


6 


5% 


8.68 


8.00 


8.44 


8.16 


7.58 


8.38 






10% 


13.46 


12.14 


12.88 


12.62 


11.76 


12.10 




10 


5% 


7.28 


6.94 


7.40 


6.42 


6.00 


7.00 






10% 


12.34 


11.50 


12.04 


11.40 


11.00 


11.14 




15 


5% 


6.86 


6.66 


6.80 


5.80 


5.44 


6.04 






10% 


11.84 


11.30 


11.50 


10.40 


10.16 


10.80 


128 


7 


5% 


7.78 


7.16 


7.58 


7.44 


6.96 


7.7 






10% 


12.30 


11.68 


11.94 


11.68 


11.22 


11.28 




12 


5% 


7.28 


7.08 


7.16 


6.52 


5.88 


6.4 






10% 


11.40 


11.12 


11.20 


10.84 


10.18 


10.46 




20 


5% 


6.82 


6.70 


6.98 


5.86 


5.50 


5.8 






10% 


11.38 


11.04 


11.26 


10.30 


9.88 


10.24 



(b) 








Tnie) 




Tni'l) 


n 


Bn 


a% 


BAR 


TUK 


PAR 


BAR 


TUK 


PAR 


64 


6 


5% 


11.22 


10.64 


9.98 


7.36 


7.02 


7.22 






10% 


15.82 


14.92 


14.40 


11.70 


10.82 


11.06 




10 


5% 


11.86 


11.66 


10.86 


7.12 


6.72 


6.94 






10% 


17.26 


17.02 


15.46 


11.68 


11.44 


11.50 




15 


5% 


12.66 


12.68 


11.78 


6.68 


6.72 


6.82 






10% 


18.36 


17.88 


17.28 


11.80 


11.28 


11.52 


128 


7 


5% 


11.40 


10.80 


10.50 


6.64 


6.48 


6.60 






10% 


16.18 


15.48 


14.56 


10.90 


10.40 


10.34 




12 


5% 


12.50 


12.00 


11.48 


6.36 


6.38 


6.20 






10% 


17.30 


17.14 


16.04 


10.34 


10.44 


10.48 




20 


5% 


13.52 


13.22 


12.48 


5.38 


5.42 


5.98 






10% 


19.24 


18.96 


17.58 


10.22 


10.26 


10.34 
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Table 2: Rejection rates in percentage under Alternative 1: (a) when the data are 
generated from model pT|) . (b) when the data are generated from model (fT2|) . 



(a) 








Tn{e) 






n 


Bn 


a% 


BAR 


TUK 


PAR 


BAR 


TUK 


PAR 


64 


6 


5% 


84.00 


81.38 


84.44 


82.72 


79.28 


83.02 






10% 


89.56 


87.22 


89.78 


88.62 


86.22 


89.00 




10 


5% 


79.04 


73.04 


79.34 


77.06 


72.30 


77.32 






10% 


85.62 


82.26 


86.08 


84.66 


80.38 


85.52 




15 


5% 


72.46 


65.80 


72.58 


70.20 


63.64 


71.80 






10% 


80.44 


75.94 


81.56 


80.00 


74.04 


79.66 


128 


7 


5% 


99.14 


98.56 


99.28 


98.90 


98.38 


98.96 






10% 


99.58 


99.36 


99.58 


99.54 


99.20 


99.56 




12 


5% 


98.00 


96.80 


98.06 


97.94 


96.60 


97.84 






10% 


99.00 


98.40 


99.16 


98.90 


98.44 


99.04 




20 


5% 


95.90 


93.12 


96.12 


95.50 


92.70 


95.70 






10% 


98.00 


96.58 


98.02 


97.60 


96.22 


97.80 



(a) 








Tn{0) 




Tni'l) 


n 


Bn 


a% 


BAR 


TUK 


PAR 


BAR 


TUK 


PAR 


64 


6 


5% 


82.06 


77.92 


83.28 


85.50 


81.66 


86.48 






10% 


89.02 


86.46 


90.40 


90.90 


88.76 


91.52 




10 


5% 


73.92 


68.04 


74.48 


78.90 


72.94 


79.24 






10% 


83.20 


79.34 


84.40 


86.16 


81.64 


86.60 




15 


5% 


65.72 


58.20 


66.54 


71.44 


63.02 


71.64 






10% 


76.70 


71.10 


78.34 


80.60 


74.80 


81.14 


128 


7 


5% 


98.84 


97.80 


99.02 


99.42 


99.10 


99.52 






10% 


99.50 


99.24 


99.66 


99.76 


99.60 


99.76 




12 


5% 


97.20 


95.20 


97.24 


98.72 


97.62 


98.86 






10% 


98.70 


97.90 


98.96 


99.44 


98.94 


99.44 




20 


5% 


93.50 


90.02 


93.74 


96.96 


94.42 


97.02 






10% 


96.98 


94.86 


97.38 


98.38 


97.12 


98.66 
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Table 3: Rejection rates in percentage under Alternative 2: (a) when the data are 
generated from model pT|) . (b) when the data are generated from model (fT2|) . 



(a) 








Tn{e) 






n 


Bn 


a% 


BAR 


TUK 


PAR 


BAR 


TUK 


PAR 


64 


6 


5% 


61.48 


63.48 


57.98 


58.28 


59.88 


55.26 






10% 


73.80 


73.98 


70.12 


70.68 


71.54 


67.90 




10 


5% 


68.04 


66.98 


65.06 


64.22 


64.08 


61.58 






10% 


78.56 


78.48 


76.42 


75.12 


74.80 


73.58 




15 


5% 


71.40 


69.62 


68.16 


66.56 


66.04 


65.20 






10% 


81.24 


80.76 


79.10 


79.00 


77.32 


75.52 


128 


7 


5% 


95.58 


95.36 


92.76 


93.60 


93.68 


91.48 






10% 


97.98 


98.08 


96.34 


97.12 


97.14 


95.38 




12 


5% 


98.52 


98.02 


96.94 


98.02 


97.32 


95.42 






10% 


99.52 


99.32 


98.54 


99.20 


98.92 


98.04 




20 


5% 


99.36 


99.30 


98.90 


98.82 


98.98 


98.38 






10% 


99.80 


99.82 


99.74 


99.64 


99.66 


99.42 



(a) 








Tn{0) 




Tni'l) 


n 


Bn 


a% 


BAR 


TUK 


PAR 


BAR 


TUK 


PAR 


64 


6 


5% 


78.10 


76.14 


76.40 


77.82 


75.46 


75.12 






10% 


86.44 


84.76 


84.54 


85.98 


84.72 


83.86 




10 


5% 


76.50 


72.84 


75.60 


77.58 


73.60 


75.92 






10% 


86.80 


84.08 


85.20 


86.78 


84.76 


85.58 




15 


5% 


75.62 


71.54 


73.44 


76.80 


72.06 


74.22 






10% 


86.98 


83.38 


84.22 


87.18 


84.00 


85.28 


128 


7 


5% 


97.76 


97.52 


100 


99.20 


98.94 


98.64 






10% 


99.46 


99.40 


99.16 


99.68 


99.60 


99.48 




12 


5% 


98.44 


97.34 


100 


99.50 


99.12 


99.04 






10% 


99.54 


99.24 


99.34 


99.80 


99.72 


99.68 




20 


5% 


98.48 


98.06 


100 


99.68 


99.58 


99.44 






10% 


99.62 


99.48 


99.40 


99.84 


99.84 


99.78 
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Table 4: Rejection rates in percentage under Alternative 3: (a) when the data are 
generated from model pT|) . (b) when the data are generated from model (fT2|) . 



(a) 








Tn{e) 




Tni'y) 


n 


Bn 


a% 


BAR 


TUK 


PAR 


BAR 


TUK 


PAR 


64 


6 


5% 


17.98 


16.64 


6.74 


13.34 


11.62 


4.28 






10% 


31.48 


27.30 


13.26 


25.8 


22.62 


9.48 




10 


5% 


37.06 


42.36 


27.24 


33.22 


40.26 


22.10 






10% 


51.16 


57.38 


41.38 


47.6 


53.58 


37.88 




15 


5% 


42.30 


46.98 


43.32 


37.90 


44.36 


41.28 






10% 


55.54 


60.26 


58.16 


53.1 


57.28 


54.12 


128 


7 


5% 


71.62 


70.66 


22.38 


65.22 


65.08 


13.44 






10% 


83.18 


83.00 


37.24 


80.44 


80.70 


28.02 




12 


5% 


85.76 


89.52 


82.76 


84.68 


89.16 


80.72 






10% 


92.78 


94.44 


91.50 


91.52 


94.10 


90.40 




20 


5% 


86.80 


88.84 


89.66 


85.00 


87.22 


89.04 






10% 


92.60 


93.50 


94.36 


91.92 


93.12 


93.80 



(a) 








Tn{0) 




Tni'y) 


n 


Bn 


a% 


BAR 


TUK 


PAR 


BAR 


TUK 


PAR 


64 


6 


5% 


17.38 


14.86 


6.46 


13.86 


11.08 


3.52 






10% 


32.18 


28.62 


14.04 


28.1 


23.16 


8.10 




10 


5% 


33.10 


38.64 


23.14 


33.44 


40.12 


22.60 






10% 


50.18 


55.18 


40.86 


49.6 


55.48 


37.98 




15 


5% 


38.20 


42.52 


38.90 


39.38 


44.38 


40.46 






10% 


52.96 


56.98 


55.48 


54.6 


60.26 


56.56 


128 


7 


5% 


63.60 


64.06 


22.40 


70.72 


69.52 


13.00 






10% 


80.98 


80.50 


39.56 


84.36 


83.54 


28.66 




12 


5% 


81.02 


85.44 


77.04 


86.74 


90.16 


84.06 






10% 


90.24 


92.64 


88.92 


93.02 


94.70 


91.56 




20 


5% 


81.54 


83.86 


85.44 


87.58 


89.08 


90.18 






10% 


89.80 


91.46 


92.44 


92.90 


93.88 


94.52 
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